Validity of Feynman's prescription of disregarding the Pauli principle 

in intermediate states 
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Regarding the Pauli principle in quantum field theory and in many-body quantum mechanics, 
Feynman advocated that Pauli's exclusion principle can be completely ignored in intermediate states 
of perturbation theory. He observed that all virtual processes (of the same order) that violate the 
Pauli principle cancel out. Feynman accordingly introduced a prescription, which is to disregard 
the Pauli principle in all intermediate processes. This ingeneous trick is of crucial importance in the 
Feynman diagram technique. We show, however, an example in which Feynman's prescription fails. 
This casts doubts on the general validity of Feynman's prescription. 
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I. INTRODUCTION 

In his space-time approach to quantum electrodynam- 
ics, Feynman advocated: "It is obviously simpler to dis- 
regard the exclusion principle completely in the interme- 
diate states" Q . He examined processes involving several 
particles and observed that all virtual processes that vi- 
olate Pauli's exclusion principle (formally) cancel out. It 
is understood that all virtual processes of the same orders 
are taken into account. On the basis of this observation 
Feynman introduced a prescription that is to disregard 
the Pauli principle in all intermediate states. This in- 
genious trick was crucial in accomplishing the enormous 
simplification and transparency of perturbation theory. 
For example, the vacuum polarization can be related to 
Feynman diagrams with an electron loop or loops. Al- 
though the process represented by a loop diagram may 
(at least partially) violate Pauli's exclusion principle, no 
restriction needs to be imposed on integrations with re- 
spect to associated momentum variables. Feynman's pre- 
scription is also often used in perturbation calculations 
for many-body systems in quantum mechanics. 

Various aspects of Feynman's prescription have been 
discussed by several authors [|j . There are some intrigu- 
ing implications regarding the meson effects in nuclei or 
nuclear matter. Feynman's prescription is instrumental 
in proving Goldstonc's theorem for many-body systems. 
We are not going to review these topics in this paper but 
we emphasize that no suspicion seems to have ever been 
raised in the literature against the validity of Feynman's 
prescription. 

The purpose of this paper we present an example that 
casts doubt about the general validity of Feynman's pre- 
scription. The example is concerned with the second or- 
der energy shift of a rclativistic bound system. We con- 



sider a model that consists of a particle bound in a given 
potential. The wave function of the particle is subject to 
the Dirac equation with the given binding potential. In 
addition to the bound particle, there is a vacuum back- 
ground. It is understood that the vacuum background 
is an integral part of the bound system. When an ex- 
ternal perturbation is applied, the energy of the system 
is shifted. We calculate the energy shift in second order 
perturbation theory. We are particularly interested in the 
vacuum effect to which the Pauli principle is relevant. 

We consider two methods, I and II, for calculating the 
energy shift. In method I we take account of the Pauli 
principle whenever it is applicable. In method II we dis- 
regard the Pauli principle altogether. We confirm that 
these two methods formally agree. This illustrates Feyn- 
man's prescription. When methods I and II are explicitly 
worked out for the example, however, the results of the 
two methods turn out to disagree with each other. We 
analyze the intriguing mechanism of this discrepancy. 

In Sec. II we set up the model and illustrate Feynman's 
prescription. In Sec. Ill we make the model more ex- 
plicit. We consider a charged particle that is bound in an 
infinite square-well potential of the Lorentz scalar type. 
This is a one-dimensional version of the "bag model" . 
For the external perturbation we assume a homogeneous 
electric field. Then the second order energy shift is re- 
lated to the electric polarizability of the system. We 
carry out the calculations of methods I and II. The two 
methods result in different energy shifts (and hence dif- 
ferent values of the electric polarizability). We analyze 
the source of the discrepancy. In Sec. IV we confirm the 
result of method II by repeating the calculation by using 
the Dalgarno- Lewis (DL) method |)-|. A summary and 
discussions are given in Sec. V. Some details concerning 
the series that appear in method II are relegated to the 
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Appendix. 



II. FEYNMAN'S PRESCRIPTION 

As a way of setting up notation, let us start with a 
problem of single-particle quantum mechanics. Let the 
Hamiltonian of the model be 

H = H + V, (1) 

where Hq is the Dirac Hamiltonian with a binding po- 
tential and V is the external perturbation. (Imagine 
something like a hydrogen atom, with Hamiltonian Hq, 
placed in a weak external electric field V. Assume that 
the proton is merely a source of the Coulomb potential 
that binds the electron of the atom.) We take Hq as the 
unperturbed Hamiltonian and treat V by perturbation 
theory. It is understood that the solutions of the Dirac 
equation with Hq are known for all stationary states, 

H \i) = ei\i), H \-j) = e- j \-j), (2) 

where i = 1, 2, ■ • • and —j = — 1, —2, • • • . The \i)'s are 
positive energy states with a > and | — j')'s are neg- 
ative energy states with e-j < 0. In particular |1) is 
the lowest positive energy state (like the lS'-state of the 
hydrogen atom). We are assuming that the eigenvalues 
are all discrete. (It is not difficult to include continuum 
states; see Sec. V.) The |i)'s and | — j)'s form a complete 
orthonormal basis set. Figure 1 schematically shows the 
unperturbed energy spectrum. 
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FIG. 1. Unperturbed energy levels defined by Eq. (Q). 
The dots indicate occupied states. 

For the unperturbed state let us take |1), the state of 
the lowest positive energy || . In single-particle quantum 
mechanics we do not consider the vacuum background. 
In other words, we assume that the negative energy states 
are all empty. In this sense let us momentarily ignore the 
dots of the states of the | —j) 's in Fig. 1. Then the energy 
shift Wqm of state |1) caused by perturbation V is given 
by 

Wqm , e ^l +e iy^l, (3) 

rrf ei - e; ^ ei - e-j 

*#1 3 



where V ia = (i\V\l) and V-j,i = (-j|V|l). Suffix QM 
refers to single-particle quantum mechanics. The sum- 
mation for i (j) is for the positive (negative) energy in- 
termediate states. 

Let us now take account of the vacuum background 
that accompanies the particle bound in |1). Let us ex- 
amine the vacuum from the point of view of Dirac's hole 
theory. The vacuum consists of an infinite number of 
particles that occupy all of the negative energy states as 
indicated with dots in Fig. 1. Before the external per- 
turbation V is applied, the negative energy states are 
eigenstates of Hq. The energy of the unperturbed sys- 
tem, including the vacuum background, is 

E = e 1+ J2^- (4) 

j 

The sum over the negative energy states is the energy 
of the vacuum. This vacuum is different from the "free 
vacuum". The ey's are eigenvalues of Hq that contains 
the binding interaction. The summation of Eq. (^) di- 
verges but this is not a serious problem. If we subtract 
the energy of the free vacuum, the summation will con- 
verge but let us not delve into this aspect. We are only 
interested in the energy shift rather than the energy it- 
self. All e's and hence the total energy E are shifted by 
external perturbation V . Let us examine two methods, I 
and II, for calculating the energy shift. 

Method I. In this method we take account of the Pauli 
principle in intermediate states. Energy e\ is shifted by 

^ = £^, (5) 

where the Pauli principle excludes the negative energy 
states as intermediate states. The difference between W\ 
and Wqm is due to the presence of the vacuum back- 
ground in the former. On the other hand the vacuum 
energy itself also shifts by 

Wvac = ^=^M ( 6) 

3 J 

Again the summation over i (j) is for positive (negative) 
energy states. The intermediate state of i = 1 is excluded 
because it is already occupied. Then the total energy 
shift of the system, including the vacuum background, is 
given by 

W = Wl + I^vac • (7) 

The Wi and W vac are both negative and hence W is 
negative. Note that Wqm can be positive because of 
the contributions from the negative energy intermediate 
states. 

Method II. In this method we disregard the Pauli prin- 
ciple in all intermediate states. We can rewrite the above 
W by adding and subtracting the same terms as 
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w' = w[ + w^ W[ = W Q 



M • 



(8) 



E 



|V- fc ,-,f 



(9) 



where | — fc)'s are also negative energy unperturbed states 
0,^). The restriction i ^ 1 has been removed in the i- 
summation for . Each of W[ and is the energy 
shift in the context of single-particle quantum mechanics. 
The W^ ac is the vacuum energy in the absence of the 
particle in The W[ and W^ ac both contain terms 
that violate the Pauli principle but such terms all cancel 
out when they are added to obtain W' . Note also that 
the effects of transitions between negative energy states 
cancel out, 



EE 

3 kjij 



|V_ 



e-fc 



0. 



(10) 



The formal equality between W of Eq. (Q) and W of Eq. 
(0) illustrates Feynman's prescription Q. Let us warn, 
however, that this equivalence relies on the convergence 
of the series involved, in particular, the j-summation of 
Wl ac of Eq. (§) that involves Eq. ©. 

In quantum field theory no negative-energy parti- 
cles appear but antiparticles of positive energies appear 
instead. The unperturbed state that we consider is 
ci'|vac). Here |vac) is the state that contains neither 
particles nor antiparticles at all. The energy of this un- 
perturbed vacuum is set to zero. The c^ is an operator 
that creates a particle with energy t\ and wave func- 
tion associated with |1). The |vac) and ci^vac) are the 
ground states of the unperturbed system within the zero- 
particle and one-particle sectors, respectively. Note that 
the particle number is a conserved quantity of the model 
under consideration. The external electric field leads to 
creation of a particle-antiparticle pair, and so on. In this 
way the whole language of hole theory can be transcribed 
into that of quantum field theory. 



III. ONE-DIMENSIONAL BAG MODEL 

We explicitly illustrate what we have shown in Sec. 
II by means of model calculations. Let us consider the 
one-dimensional bag model H[To||, which is a relativistic 
version of the infinite square- well potential model of non- 
rclativistic quantum mechanics. We define the model by 
the Dirac equation in one dimension, 

H tl)(x) = [ap + f3m + /3S(x)]<if>(x) = etp(x) , (11) 

where to is the mass of the particle, S(x) is a Lorentz 
scalar potential and p — —id/dx. We use units such 
that c = h = 1. For the 2x2 Dirac matrices, we use 



a = <7 y and /3 = cr z , where o~ y and o~ z are the usual Pauli 
matrices. For S(x), we assume that 



S(x) 



So 




for 
for 



\x\ > a 
\x\ < a 



(12) 



where So is a positive constant. It is understood that we 
let So — ► oo. For the bag model in three dimensions as a 
model of hadrons, see Ref. p|JIo|. 

With the specific choice of a, no complex numbers ap- 
pear in the Dirac equation. We write the wave function 
ip{x) as 



ip(x) 



u{x) 
v{x) 



(13) 



The u{x) and v(x) vanish outside the bag, i.e., for \x\ > a, 
and are discontinuous at |a;| = a. They are subject to the 
boundary condition 



u(±a) = ^fv(±a) 



(14) 



The scalar density i(j>{3ip = u 2 — v 2 vanishes at |a:| = a, 
but the vector density ip'tp = u 2 + v 2 does not have to 
vanish at \x\ — a. 

The solutions of Eq. ([n]) can be classified in terms of 
parity. For even parity, we obtain 



u(x) — N cos kx , v(x) 



k sin kx 

-N (15) 

e + to 



where k — ve 2 - to 2 and iV is a normalization factor. 
For negative parity, we similarly obtain 



u(x) = N sin kx , v(x) = N 



k cos kx 



Equation (|l4|) leads to 



tan ka = 



±e + m 



(16) 



(17) 



Where the double sign is + (— ) for positive (nega- 
tive) parity. Equation (^) determines k n and e„ (with 
n = 0, 1, 2, ...) for each parity. The e n can be positive or 
negative. The normalization factor N is given by 



_> _ e(e + to) 



N z = 



2ae 2 ' 



(18) 



which applies to both of Eqs. ([l5|) and (|l6|). When the 
potential for the Dirac equation is a pure Lorentz scalar, 
there is symmetry between positive and negative ener- 
gies. This symmetry is manifest in our model. For a 
positive parity state with energy e, there exists a nega- 
tive parity state of energy — e. This can be seen through 
Eq. @. 

The special case of to = is very simple and instruc- 
tive. In this case the solutions k n (with n = 0, 1, 2, ...) of 
Eq. © are given by 
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for 



even parity, 
odd parity, 



e > 0, 
e < 0, 



W 



±PS> 



(22) 



* ) 1 for 
47 a 



even parity, e < . 
odd parity, e > . 



In Sec. II we designated the energy levels with i, —j 
and —k. For the bag model, however, we denote the 
levels with n ps where p stands for parity, s is the sign 
of the energy and n = 0, 1, 2, .... For example ++ and 
h are the lowest and the second lowest positive energy 
states, respectively. They are |1) and |2), respectively, in 
the notation of Sec. II. If we denote the energy of state 
n ps with e(n ps ) we obtain 



(19) The coefficient P defines the electric polarizability of the 
system. The W can be written similarly with P' . 

Interaction V(x) connects states of opposite parity. 

The matrix element between states n +s and n' s is given 
by 



(n'-'\V\n+')=V n> . 



XNN' 



cos kx sin k x — 



kk' sin kx cos k'x 



(e + m)(e' + m) 



xdx , (23) 



e(n++) 
e(n+-) 




e(n 




where k and k! are associated with the states of n and 
n', respectively, and similarly for N, N', e and e'. Let 
us first examine the simple case of m = 0, which gives us 
much insight into the problem. Then the matrix clement 
becomes 



(20) 



\{n'- s \V\n + *)\ = 



A 



a(e-e') 2 



(24) 



Figure 2 shows the energy spectrum of this case of m = 0. 
The energy levels are all equally spaced. If mass m is 
increased from 0, the energy levels are pushed away from 
E = 0, the levels nearer to E = being more affected 
than those further away from E = 0. 



For the quantities of Eqs. (|^-0) of method I we obtain 

w = w a++ + w mc , 



W 0++ = J2f(n' 



Even Parity 



Odd Parity 



(2-+ ^) 



E = 



(2++,|) 



FIG. 2. Energy levels of the bag model with m = 0. The 
pair of numbers in brackets are respectively quantum number 
n ps and e in units of n/a. The dots indicate occupied levels. 

Let us now assume the external perturbation 

V(x) = Xx, (21) 

and work out the second order energy shifts W and W' 
explicitly for the bag model. If the charge of the particle 
is q and the external perturbation is due to a constant 
electric field £ along the x-axis, then A = — q£. Then the 
second order energy shift W takes the form 



W 



E 

71 = 



W n+ - = J2 f{n + n' + -), 



W„ — 



n'=0 



£/(„ + „'+-) 

n'=l 



W n+ - 



where 



m = 



A 2 a 3 1 



(25) 
(26) 

(27) 

(28) 

(29) 

(30) 



\ is due to 



In Eq. (26) the term with the argument n' 
the transition ++ — > n' h . In Eq. ( p^ ) the term with 
n + n' + | is due to — ► n' + . Note that n' = is 
excluded in the summation for W„.— . The above series 
all converge very rapidly. In Eq. (|26| ) the first term with 
n' = constitutes 99% of the sum. 

Next, let us turn to the W of method II, Eq. (§). Let 
us again consider the m = case. Curiously enough, it 
turns out that the energy shifts of the individual levels 
all vanish and consequently the total energy shift is zero 
in this case, i.e., in the notation of Sec. II, 
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W[ = W'_ 3 ; = , W' = . 



For example, we find that Wg ++ (= W[) is of the struc- 
ture 



K ++ = E 



n'=0 



= 0. 



(32) 



Recall that f(x) is an odd function. The two terms in the 
square brackets respectively correspond to the two terms 
of the right hand side of W QU (= W{) of Eq. (§). The 
first (second) term is due to the intermediate states of 
positive (negative) energies. Exactly the same situation 
is found for the energy shift of each of other states, that 
is, 



<+- = E 



n'=0 



/(„ + n'+|)+/(n-n' + i; 



(31) W into W; see Eq. @. In the way as IF' is explicitly 
worked out as shown above in method II, that is, the n! 
summation is done before the n summation, the series of 
Eq. (35) is actually arranged such that its sum assumes 
a nonzero value. This is where the discrepancy between 
methods I and II stems from. Let us emphasize that the 
conspiracy of the above alternating series is well hidden 
in the sense that all the (single) series that appear in the 
steps of method II are absolutely convergent. 

We have also examined the case with nonzero values of 
mass m. The calculation is lengthy but straightforward. 
So we do not describe it. We have confirmed that essen- 
tially the same situation persists, that is, the results of 
the two methods disagree. For the states of very large 
values of j and/or k, effects of the finite mass m becomes 
negligible. Therefore the nonconvergence aspect of the 
( I series involved is not essentially affected by m. Figure 3 

shows W (solid line) and W (dashed line) as functions 
of ma. Note that the difference between the two is larger 
(33) for smaller ma. Figure 3 also shows the nonrelativistic 
limit (dotted line) that we derive in Sec. IV. 



K- = E 



n'=0 



f{n + n' + \)+f{n-n'-\) 



= 0. 



(34) 



The vanishing of W' + _ and W' n __ given above may not 
be immediately obvious. In Appendix we show that the 
series of Eqs. (33) and ( |34| ) can be rewritten exactly 
in the form of Eq. @. The result of W = is in 
clear contradiction with W of Eq. ( p5|) that is nonzero 
and negative. This is very puzzling. The vanishing of 
W' means that the system is rigid against the external 
perturbation, which we find intuitively strange. 

In Sec. II we warned that the equivalence between W 
and W 1 relies on the assumption that the summations 
involved converge. What happens in the above puzzle 
is the following. Let us first explain it by using the no- 
tation of Sec. II. There is no problem in convergence 
of the summations except for the j summation of Eq. 
(^) . Each of W_j is well defined, but when it is summed 
with respect to j, we obtain the double sum of the left 
hand side of Eq. (flO|). As can be seen from Eq. (24), 
|V"-fc,-j| 2 depends on j and k only through the difference 
j — k. Coming back to the notation of this section, the 
left hand side of Eq. (|l0|) becomes 



EE 

n=0n'=0 L 



(35) 
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FIG. 3. Relativistic energy shifts W of method I (solid 
line), W' of method II (dashed line) and their nonrelativistic 
counterpart 6nr' 2 ' of Eq. (^) (dotted line) are shown in 
units of A 2 a 3 as functions of ma. The £nr' 2 ' is not meaningful 
unless ma 3> 1. 



IV. THE DALGARNO-LEWIS METHOD 



The summations with respect to n and n' individually 
converge. When the two summations are combined, how- 



ever, we realize that Eq. (35) involves something like the 
alternating series 1 — 1 + 1 — 1 + 1... . We show this 
explicitly in Appendix. This series can converge only 
conditionally at best. Its sum depends on how the series 
is arranged. Equation (pq) was set to zero in rewriting 



The calculation of method II that was presented in 
Sees. II and III is somewhat involved. So it would be 
good to confirm it by repeating the calculation in a dif- 
ferent manner. We do so by using the DL method. The 
DL method is an alternative form of perturbation the- 
ory in which summations over intermediate states are 
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avoided [|3|-|5|]. As a price for it, one has to solve an inho- 
mogeneous differential equation. The DL method is often 
used in calculating the electric polarizability of nonrela- 
tivistic bound systems. There is another similar, pow- 
erful method called logarithmic perturbation expansion 
0,^2|, which we do not use here. Consider any one of 
the energy levels. Let its unperturbed wave functions be 
^°\x) and its first order perturbation be ip^ix). The 
^^(x) can be determined by the DL equation 



(Hq - e (0) )^>« (x) = -V(x)iP {0) (x) , 



(36) 



where e^ ^ is the unperturbed energy, i.e., e of Sec. III. 
We write ip^ as 



ib (r Ux) - ( uW S x) 

* [X) ~[vW(x) 



(37) 



Its components are again subject to boundary condition 
(0), 



v (1 \a) = -«W( ) 



(38) 



When V^ 1 ) is found, the second order energy shift e^ 2 ' 
can be calculated by 



.(2) 



V{x)ifj ni) (x)^°\x)dx. 



(39) 



The summation over intermediate states is done implic- 
itly. It is understood that there is no restriction on in- 
termediate states due to the Pauli principle. Since it 
includes all intermediate states of negative as well as 
positive energies, is nothing but W[ or W'_j of Sec. 
II. If we take the unperturbed wave function for |1) for 
\ip^(x), for example, we obtain W[. Obviously the DL 
method is not useful for method I. 
By solving Eq. (pq) we obtain 



(x) 



v {1) (x) = 



XN 
2k 2 



XmN 
2ek 2 



[mi cos kx + r\tk(x 2 — a 2 ) sin kx 



rjmk 
~2~T 



kx sin kx + 

(x 2 - a 2 ) 



2a sin kx 



(40) 



e + m 

e 2 (e — m) 



(e — m)a 

2 k ; 



m 



cos ka 



(41) 



where 77 = 1( — 1) for even (odd) parity. The e and k are 
e^ ) and k^ , respectively. With these and in 
Eq. ( |39| ) we arrive at 

(2) _ A^m 



24k 4 ea[2{ea) 2 + ma] 
x {2(ka) 2 (ma + 3)[4(ea) 2 - 6ma - 3] 

~15(ma) 2 (2ma+ 1)} . (42) 



This applies to any of the energy levels with an appro- 
priate choice of k that is subject to Eq. (0). We have 
explicitly confirmed that agrees with W{ or W'_j of 
method II for each of the energy levels and hence the 
same W as that of Sec. III. The vanishes if m = 0. 
This is consistent with what we found in Sec. III. 

In the non-relativistic limit of m — » 00, Eq. d38| ) is 
reduced to 



CNR 



(2) 



Am 
24k 4 



[4(ka) 2 - 15] 



ka = — 
2 



(43) 



The nonrelativistic value of ka follows from Eq. ( |17j ) 
with the plus sign and m — > 00. Equation (|3^) agrees 
with the result for an infinite square-well potential of 
non-relativistic quantum mechanics M]. The eNR^ is 
compared with its relativistic counterparts W and W' in 
Fig. 3. Note that even when m is as large as m = 3/a (or 
m w 600 MeV if a = 1 fm), the relativistic energy shifts 
are about twice as large as their nonrelativistic counter- 
part. 



V. SUMMARY AND DISCUSSIONS 

For a system consisting of a particle bound in a given 
potential together with its vacuum background, we ex- 
amined the second order energy shift caused by external 
perturbation V . We examined two formally equivalent 
methods of calculation, I and II. Method I takes account 
of the Pauli principle in intermediate states whenever it 
is applicable. In method II the Pauli principle is com- 
pletely ignored. We showed that, if the energy shifts 
of all occupied levels are summed up in method II, the 
terms violating the Pauli principle formally cancel out. 
Thus the two methods appear equivalent. This illustrates 
Feynman's prescription. 

This equivalence, however, is not free from ambigu- 
ity. We calculated the energy shift explicitly for the 
one-dimensional bag model with external perturbation 
V(x) — Xx. As shown in Fig. 3, the two methods lead 
to different energy shifts. Thus Feynman's prescription 
fails in this example. For method II, we did the calcula- 
tion in two different manners, one by summing up over 
the intermediate states and the other by using the DL 
method. The same result were obtained by two calcula- 
tions. In method II the energy shifts of the individual 
occupied levels are unambiguously obtained. When they 
are summed over all negative energy states, however, an 
alternate series emerges. The sum of the series depends 
on how the summation is done. This is essentially the 
source of the discrepancy between the two apparently 
equivalent methods. The alternating series is hidden such 
that, if one simply follows method II, one would not no- 
tice it. 

Feynman's prescription fails in the specific example 
that we have described. A question naturally arises here. 
Does similar difficulty arise in more general situations? 
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We suspect that it may well do. Let us first point out 
that, although we assumed a specific form of external 
perturbation V(x) = Xx, Feynman's prescription fails in 
the one-dimensional bag model irrespectively of the form 
of V(x). Again for simplicity let us assume m = 0. Then 
the matrix element of V{x) are of the form of 

/a />a 
V(x) sin[(fc — k')x]dx or / V(x) cos[(fc — k')x]dx . 
-a J —a 

(44) 

No matter how large k and k 1 become, the matrix element 
is of the same form as that of Eq. (p4]). This feature re- 
mains essentially the same when m becomes nonzero. In 
this connection, recall what we said in the last paragraph 
of Sec. III. 

Next let us consider the three-dimensional bag model, 
subject to a constant external electric field. The pertur- 
bation interaction can be taken as Xz = XrcosO. For 
states with large quantum numbers, the radial part of 
the wave function is similar to the one-dimensional wave 
function. This is so in the sense that at large distances 
the spherical Bessel functions involved are like the sine 
and cosine functions. For the angular part, the matrix el- 
ement of cos (9 between two adjacent angular momentum 
states has a part that remains finite no matter how large 
the angular momenta become. Therefore, the alternating 
series involved in method II will remain. We are aware 
of a few calculations of the electric and magnetic polar- 
izabilities of the nucleon by using the bag model Q |l6) . 
Method I was used in these calculations and hence the 
problem with method II was not encountered. 

We have assumed that the energy spectrum of the un- 
perturbed system is discrete. The case of continuum 
spectrum can be handled by enclosing the system in a 
very large cavity. The unperturbed Hamiltonian Hq in 
this case can be that of the bag model (with a large ra- 
dius) plus some other interaction that produces states 
localized, say, around the origin. Let us consider such a 
model in one dimension. The perturbation of the form 
of V(x) = Xx, if taken literally, would not make much 
sense because it becomes very large as x approaches the 
cavity radius. If one chooses V(x) such that it remains 
reasonably small within the entire cavity, one can treat 
it by perturbation theory. Then the calculation will go 
in a way essential the same as we have done. Feynman's 
prescription will probably fail again. 

As far as we know the example that we have presented 
is the first counter-example against Feynman's prescrip- 
tion that seems to have been taken for granted for many 
years. If we have to choose between methods I and II, 
we are inclined to take method I and abandon method II 
that is based on Feynman's prescription. We think that, 
if we encounter ambiguity by disregarding the Pauli prin- 
ciple, we should remain faithful to the Pauli principle in 
every step of calculation. In view of the fact that Feyn- 
man's prescription has been used extensively, its possible 
failure may have serious implications. 
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Appendix 

Let us first examine how W + _ of Eq. ( |3^ ) and W' of 
Eq. ( |34| ) vanish. First note that 

K+- = W[ n+1) - ■ (45) 

Therefore it is sufficient to show that W = 0. This 

can be seen as follows: 

W' n - = £ f(n + n> + 1) + 2 /(n - n' - ±) 

ti'^O n'=Q 
oo 1 

+ J2f(n-n'--). (46) 

n'— n 

If we define v — n' — n, the last sum can be reduced to 

oo 1 

E/(-^-o)- ( 4? ) 

i/=0 

It is not difficult to see that the first two sums can be 
combined into 

oo 1 

E/("+ 5 )- ( 48 ) 



Because of /(x) = —f(—x), W' n __ = follows. Al- 
though W'— can be regarded as an alternating series, 
it is absolutely convergent. It is not like the alternating 
series that we mention below Eq. (^5|). 

Next let us examine the structure of the double series 
ofEq. ((3^). Consider a set of (n, nf) such that n = n'+l, 
i.e., 

(n,n') = (1,0), (2,1), (3,2), ••• . (49) 

For this set we find that the term in the square brackets 
of Eq. (H|) takes the same value /(§)+/(!)• ^ ms ^ s so 
no matter how large n and n' individually are. Similarly, 
for a set of (n, n') such that n = n' — 1, i.e., 

(n,n') = (0,l), (1,2), (2,3), ••• , (50) 

we find /(-i)+/(-§) = -[/(§) + /(±)]. Therefore, the 
terms corresponding to the sets of (n,n' = n ± 1) can be 
seen as an alternating series like 1 — 1 + 1 — 1 + 1 ■ ■ ■ . If 
we pair the above like (1,0) and (0,1), (2,1) and (1,2), 

• • •, then we find the double sum vanishes, like Eq. (|lC|). 
If we pair the above like (1,0) and (1, 2), (2, 1) and (2,3), 

• • •, then the sum does not vanish. We find similar series 
for (n,n = n±2), (n, n — n±3), and so on. This shows 
that the sum of the double series has an ambiguity that 
is related to how the n-n' summation is done. 
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